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ABSTRACT 

We describe the discovery of a massive transiting hot Jupiter with a very short 
orbital period (1.30619 d), which we name TrES-3. From spectroscopy of the 
host star GSC 03089-00929, we measure T^s = 5720 ± 150 K, log^ = 4.6 ± 0.3, 
and vsini < 2kms^^, and derive a stellar mass of 0.90 ± 0.15 M©. We estimate 
a planetary mass of 1.92 ±0.23Mjup, based on the sinusoidal variation of our 
high-precision radial velocity measurements. This variation has a period and 
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phase consistent with our transit photometry. Our spectra show no evidence of 
hne bisector variations that would indicate a blended eclipsing binary star. From 
detailed modeling of our B and z photometry of the 2.5%-deep transits, we deter- 
mine a stellar radius 0.802 ± 0.046 Rq and a planetary radius 1.295 ± 0.081 -Rjup- 
TrES-3 has one of the shortest orbital periods of the known transiting exoplan- 
ets, facilitating studies of orbital decay and mass loss due to evaporation, and 
making it an excellent target for future studies of infrared emission and reflected 
starlight. 

Subject headings: planetary systems — stars: individual (GSC 03089-00929) — 
techniques: photometric — techniques: radial velocities 



Introduction 



The OGLE-III deep-field survey ( lUdalski et al.ll2003l . and references therein) has iden- 
tified 3 transiting planeta ry systems with very short orb ital periods: OGLE-TR-56b (P 



29hr: iKonacki et all 120031 ) . 0GLE-TR-11 3b {P = 34 hr 



2004h . and OGLE-TR-132b (P = 41 hr; iBouchv et al. 



Bouchv et al. 2004; Konacki et al. 



20041 ). At the time of discovery. 



such 1-d period planets were conspicuously absent from radial velocity (RV) surveys. These 
surveys had not discovered a planet with an orbital period of substantially less than 3d, 
despite the sensitivity of such surveys to planets with short orbital periods. Several au- 
thors reconciled this a pparent discrepancy by determining the respective observational biases 



and selection effects ( IGaudi et al.l l2005l : iPont et al.l l2005l : iGould et al.l l2006l : iFressin et al. 
20071 ). Since then, hot J upiters with periods just over 2d have been detected by both 



RV surveys (HP 189733b; 



Bouchy et al.l l2005l ) and wide-field transit surveys (WASP-2b; 



Collier Cameron et al.l 120071 ). indicating a consistent picture for the underlying distribution 



of orbital periods of hot Jupiters. 



With the discovery of these nearby planets, we continue to observe (ISozzetti et al.ll2007l ) 
the decreasing li near relation b etwee n P a nd the plane t ary rn ass (Mp) for transiting planets 
first noticed by iMazeh et al.l (120051 ) and I G audi et al.l (120051 ). It appears from this small 
number of planets that shorter period exoplanets are on average more massive. This may 
suggest different ev olution of these planetary systems than the canonical migration (see, e.g.. 
Trilling et al.lll998l ) of hot Jupiters from beyond the ice line. 



We present here the discovery of TrES-3, a massive planet which has one of the shortest 
periods of the transiting planets, and which was independently identified by two transit 
surveys. 
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2. Observations and Analysis 



Transits of the parent star GSC 03089-00929 were detected by two 10-cm telescopes, 
part of the Trans- at lantic Exoplanet Survey (TrES) network. Sleuth (Palomar Observatory, 
California) observed 7 full transits and 3 partial transits, of which 3 full and 2 partial 
transits were observe d by the Planet Search Survey Telescope (Lowell Observatory, Arizona; 
Dunham et al.ll2004l ). We obtained these observations between UT 2006 May 6 and 2006 
August 2, during which time we monitored a 5?7 x 5?7 field of view (FOV) centered on 
6 Her. We then processed the data from each telescope sepa rately, as we have described 
previously (IDunham et al.ll2004j : lO'Donovan et al.ll2006bl . 120071). We then searc h our binned 
light curves using the box-fitting transit search algorithm of iKovacs et al.l (120021 ) for periodic 
event consistent with the passage of a Jupiter-sized planet across a solar-like star. We flagged 
TrES-3 as a candidate, noting that the transit duration of 1.3 hr implied a high impact 
parameter h. 

Transits of TrES-3 were independently observed by the Hungarian Automated Tele- 
scope Network (HATNet; ISakos et al.ll2004l ). We observed the HAT fleld "G196" between 
UT 2005 June 8 to December 5, using the HAT-7 telescope at the Fred L. Whipple Ob- 
servatory (FLWO) and HAT-9 at the Submillimeter Array atop Mauna Kea. We applied 
st andard calibration a nd aperture photometry procedures to the frame s as described earlie r 



m 



Bakos et al.l (120071 ). We also applied the trend filtering algorit hm (iKovacs et al.ll2005l) . 



and identified TrES-3 as a transit candidate using the algorithm of iKovacs et al 



J2002h 



We observed TrES-3 using the CfA Digital Speedometers (jLathamlll992l ) on 13 occa- 
sions from 2006 September 9 to 2007 April 6. These spectra are centered at 5187 A and 
cover 45 A with a resolving power of A/AA ~ 35,000. We derived RVs at each epoch by 
cross-correlation against synthetic spectra created by J. Morse using Kurucz model stellar 
atmospheres (J. Morse & R. L. Kurucz 2004, private communication). There is no signif- 
icant variation in these measurements, which have a mean of -1-9.58 km s~^ and a scatter 
of 0.73kms~^. This limits the mass of the companion to be less than about 7Mjup. We 
estimated the stellar parameters from a cross-correlation analysis similar to that above, 
assuming a solar metallicity. The effective temperature (Teg), surface gravity (log (7), and 
projected rotational velocity (f sini) we derive are listed in Tabled! Based on these values, 
we estimate the stellar mass to be = 0.90 ± 0.15 Mq. Furthe r spectroscopic analysis o f 
the parent star (similar to that performed for TrES-1 and TrES-2; ISozzetti et al.ll2004l . 120071 ) 
is in progress and will be presented elsewhere. 

The 2MASS J - K color (0.407mag) and UCAC2 proper motion (39.1 mas yr'^) of 
GSC 03089-00929 are as expected for a nearby G-dwarf. We obtained absolute BVRcIc 
photometry for TrES-3 on the night of UT 2007 April 14 with the 105-cm Hall telescope 
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at Lowe l l Obs ervatory. The photometry was calibrated using seven standard fields from 
Landoltl (Il992l ). The results are listed in Table [H 



We observed a full t ransit in the z-band on UT 2007 March 26 using KeplerCam (see, 



e.g., iHolman et al.l l2006l ) at the FLWO 1.2-m telescope. We obtained 150 90-s exposures 
of a 23^1 X 23fl FOV containing TrES-3. We observed another transit in Bessell B on UT 
2007 April 8 with the Las Cumbres Observatory Global Telescope 2-m Faulkes Telescope 
North on Haleakala, Maui, Hawaii. The CCD camera imaged a 4.6' square field with an 
effective pixel size of 0'.'28 (the images were binned 2x2). We obtained 126 images with an 
exposure time of 60 s and a cadence of 70 s. We defocused the telescope slightly (to an image 
diameter of about 3") to avoid saturation of the brightest comparison star. We analyzed 
these data sets separately. Using standard IRAF procedures, we corrected the images for 
bias, dark current, and flat-field response. We determined fluxes of the target star and 
numerous comparison stars using synthetic aperture photometry. To correct the target star 
fluxes for time-varying atmospheric extinction, we divided them by a weighted average of 
the fluxes of all the comparison stars. 

We collected high-precision RV measurements on UT 2007 March 27-29 using HIRES 
( IVogt et al.lll994l ) and its I2 absorption cell on the Keck I telescope. We obtained a total 
of 7 star+iodine exposures and 1 template spectrum. All spectra were gathered with a 
nominal resolving power A/AA ~ 55 000, using the HIRES setup with the 0.86" slit. We 
used an integration time of 15 min, which yielded a typical signal-to-noise ratio of 110 pixel"^. 
We reduced the raw spectra using the MAKEE software written by T. Barlow. We have 
described the software used to derive relative radial velocities with an iodine cell in earlier 



works (e.g., lO'Donovan et al.ll2006al : ISozzetti et al.l 120061 ). We estimate our internal errors 
from the scatter about the mean for each spectral order divided by the square root of the 
number of orders containing I2 lines, and find them to be around 10 ms"^. The radial- velocity 
measurements are listed in Table [2J 

We constrained the orbital fit to these data to have zero eccentricity (e), as expected 
from theoretical arguments for such a short period planet, and we also held P and the transit 
epoch Tc fixed at their values determined from the photometric data. The rms residual 
from this fit (15.4ms~^) is larger than the typical internal errors (10ms~^). Preliminary 
analysis of the template spectrum suggests that the star shows evidence of activity. For the 
i nferred spectral t ype, the presence of radial- veloci ty "jitter" of 10-20 ms~^ is not unexpected 
(ISaar et al.lll998l : ISantos et al.ll2000l : IWrightl 120051 ) . and would explain the excess scatter we 
find. Figure [T] displays the RV data overplotted with the best-fit model, along with the 
residuals from the fit. The parameters of the orbital solution are listed in Table [31 We 
find a minimum planetary mass of Mp sin i = 2.035 ± 0.090[(Mp + M^)/Mq]2/3Mjup, where 
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i represents the orbital inclination. As a consistency check we repeated the fit with P fixed 
and e = as before, but solving for T^. The result is = HJD 2,454,185.911 ± 0.045, which 
is consistent with, but less precise than, the value determined from the photometry (Table [3]). 

We investigated the possibility that the velocity variations are due not to a planetary 
compan ion, but instead to distortions in the line profiles caused by an unresolved echpsing 
binary (ISantos et al.ll2002l : [Torres et al.ll2005l ). We cross-correlated each spectrum against a 
synthetic template matching the properties of the star, and averaged the correlation functions 
over all orders blueward of the region affected by the iodine lines. From this representation 
of the average spectral line profile we computed the mean bisectors, and as a measure of 
the line asymmetry we calculated the "bisector spans" as the velocity difference between 
points selected near the top and bottom of the mean bisectors flTorres et al.ll2005l ). If the 
velocity variations were the result of a stellar blend, we would expect the bisector spans to 
vary in phase with the photometric period with an amplitude similar to that seen in the 
RVs (iQueloz et al.ll200ll : iMandushev et al.ll2005l ). Instead, we do not detect any variation 
exceeding the measurement uncertainties (see Figured]). We conclude that the RV variations 
are real, and the star is indeed orbited by a Jovian planet. 



3. Estimates of Planet Parameters and Conclusions 



We analyze our photometric time series using the analytic transit curves of lMandel fc Agol 



(2002) and the Markov Cha i n Mo nte C arlo (MCMC ) techn iques described in iHolman et al. 



(120061 ) , ICharbonneau et al.l (120071 ) , and IWinn et al.l (120071 ) . We assume a circular orbit of 
constant P. We first estimate by fitting a model light curve (as described below) to only 
the z data. We then determine P by phase-folding the z data with the TrES and HAT 
discovery data (which affords a baseline of 1.8 years) while varying the trial values for P. We 
then fix the values for Tc and P (stated in Table [3]) in the subsequent analysis. 

The values of the planet radius Rp and the stellar radius i?^ as constrained by the 
light curves are covariant with M^. In our MCMC analysis, we fix = 0.9 M©, and then 
estimate the systematic error in the radii using the scahng relations Rp oc R^, oc Mj (see 
footnote to Tabled]). We assume a q uadratic limb- darken ing law with coefficients fixed at the 
band-dependent values tabulated in IClaretl (120001 . 120041 ) for the spectroscopically-estimated 
Teflf and log g and assuming solar metallicity. 

The remaining free parameters in our model are Rp, R^, and i. We require two additional 
parameters, fc^ and k^, the respective residual color-dependent extinction to the B and z 
light curves, assuming that the observed fiux is proportional to exp {—km), where m denotes 
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the airmass. We find that the TrES and HAT discovery data are too noisy to meaningfully 
constrain the parameters, and so we restrict our analysis to the B and z data. We first 
f ind the values of Rp, R^, i, ks, and that minimize using the AMOEBA algorithm 



( iPress et al.lll992l ). This model is shown as the solid curves in Figured We then created 
2 MCMC chains with 428,000 points each, one starting from the best-fit values and one 
starting a randomly-generated perturbation to these values. We subsequently rejected the 
first 100,000 points to minimize the effect of the initial conditions, and found the results of 
the two chains to be indistinguishable. We then examined the histograms of the 5 input 
parameter values, as well as the histograms for several combinations of parameters relevant 
to anticipated follow-up studies. We assigned the optimal value to be the median, and 
the l-cT error to be the symmetric range about the median that encompassed 68.3% of the 
values. We state our estimates of the parameters in Table [31 The estimated values for 
ke (0.0029 ± 0.0006) and k, (-0.0021 ± 0.0004) are small and of opposite sign, which is 
consistent with a modest difference between the average color of the calibration field stars 
and the target. 

Despite the V-shaped transit, our best-fit values indicate that TrES-3 is not grazing, 
i.e. the disk of the planet is entirely contained within the disk of the star at mid-transit, 
although grazing solutions are permitted by the data. Importantly, our abihty to obtain 
well-constrained estimates of Rp and i?* despite the large b hinged on having observations 
of the transit in both B and z. The large difference in central wavelength and hence stellar 
limb- darkening between these two bands permitted us to rule out a family of degenerate 
solutions that is allowed by observations in only a single color. We tested this notion by 
repeating the analysis above for only the z data, and found that values of Rp as large as 
2.0-Rjup could not be excluded. 

TrES-3 presents a useful testbed for theoretical models of gas giants. Its radius places 
it in the growing family of planets with radii that exceed that predicted for models of irradi- 
ated gas giants. It is one of the most massive transiting planets, and has one of the shortest 
periods. We recall that the discovery of OGLE-TR-56b at a distance of only 0.023 AU from 
its star stimulated investigations into the timescales for orbital decay and thermal evapo- 
ration. The comparable orbital separation of TrES-3 implies that many of these estimates 
are directly applicable to the new planet, but with the important difference that the much 
brighter apparent magnitude affords the opportunity for more precise study. In particular, 
we note that direct searches for decay o f the orbital pe riod may inform our understanding 



of dissipation in stellar convective zones (ISasselovl 120031 ). particularly since both the TrES-3 



planet and its stellar convective zone are more massive than that of OGLE-TR-56b. Fur- 
thermore, the mass and orbital separati on of TrES-3 are intriguingly close to the critical 



values estimated by iBaraffe et al.l (120041 ) below which evaporation would become a domi- 
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nant process. The measurement of the angle between the planetary orbit and stellar spin 
axis of TrES-3 may detect the substantial misalignment t hat might be expected for planets 



that were tidally captured rather than migrating inwards (iGaudi &: WinnI 120071 ) . Assuming 
isotropic emission, the equilibrium temperature of TrES-3 is 1643(1 — Ay/'^ K, where A is 
the Bond albedo. We intend to obtain Spitzer observations of TrES-3, as we have previously 
done for TrES-1 and TrES-2. Fina lly, we note that TrES-3 is extremely favo r able for at 



tempts to detect reflected starlight (ICharbonneau et al.lll999l : iLeigh et al.ll2003l : iRowe et al. 



20061 ) and thus determine the geometric albedo, p\ of the planet. The flux of the planet near 
opposition relative to that of the star is px x [Rp/af = Px x 7.5 x 10~^, which is more than 
twice that of any of the other known nearby transiting planets. 
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Table 1. TrES-3 Parent Star 





Parameter 


Value 


Reference 


R.A. 


(J2000) 


17 52™07!03 




Decl. 


(J2000) 


+37 32 46.1 




GSC 




03089-00929 






(Mq) 


0.90 ±0.15 


1 




(Re) 


0.802 ± 0.046 


1 




(K) 


5720 ± 150 


1 


\ogg 


(dex) 


4.6 ±0.3 


1 


vsini 


(kms"-'-) 


< 2 


1 


V 


(mag) 


12.402 ± 0.006 


1 


B-V 


(mag) 


0.712 ± 0.009 


1 


V-Rc 


(mag) 


0.417 ±0.010 


1 


V-Ic 


(mag) 


0.799 ±0.010 


1 


J 


(mag) 


11.015 ±0.022 


2 


J-H 


(mag) 


0.360 ±0.030 


2 




(mag) 


0.407 ±0.028 


2 




(mas yr~^) 


[-22.5,32.0] 


3 



References. — (1) This work; (2) 2MASS Catalog; (3) 
UCAC2 Bright Star Supplement 

^The uncertainty in R^, includes both the statistical error 
and the 5.6% uncertainty resulting from uncertainty in M^. 
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Table 2. Relative radial- velocity measurements of TrES-3 



Observation Epoch 


Radial Velocity 




HJD - 2,400,000 


m 


m s~^ 


54186.99525 


229.0 


10.8 


54187.11354 


82.1 


8.5 


54187.96330 


167.5 


8.4 


54188.04216 


286.3 


9.7 


54188.96503 


-336.5 


8.6 


54189.04672 


-234.5 


10.9 


54189.09622 


-183.9 


9.3 
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Table 3. TrES-3 Planet 



Parameter 


Value 


P (d) 


1.30619 ±0.00001 


n (HJD) 


2454185.9101 ± 0.0003 


a/R^ 


6.06 ±0.10 


a (AU) 


0.0226 ±0.0013 


b = a cos i/R^, 


0.8277 ±0.0097 


n 


82.15 ±0.21 


K (ms-i) 


378.4 ±9.9 


7 (ms~^) 


-70.6 ±6.2 


Mp (Mj,p) 


1.92 ±0.23 




0.1660 ±0.0024 


Rp ^ (-Rjup ) 


1.295 ±0.081 



^The uncertainty in Rp includes both the sta- 
tistical error and the 5.6% uncertainty result- 
ing from uncertainty in (Tabled]). 

^-Rjup = 71,492 km, the equatorial radius of 
Jupiter at 1 bar. 
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Fig. 1. — Top: Radial- velocity observations of TrES-3 obtained with Keck/HIRES using the 
I2 cell, shown relative to the center of mass and adopting the ephemeris in Table [31 The 
best-fit orbit {solid line) is overp lotted. Middle: Residuals from the best-fit model to the 
radial velocities. Bottom: Bisector spans shifted to a median of zero, for each of the iodine 
exposures as well as for the template (which is shown as the additional data point at phase 
0.937). 
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Fig. 2. — Relative flux of the TrES-3 system as a function of time from the center of transit, 
adopting the ephemeris in Table [3l and including the residual color-dependent extinction 
correction (§3]). Each of these follow-up light curves is labeled with the telescope and filter 
employed. We have overplotted the simultaneous best-fit solution, adopting the appropriate 
quadratic limb-darkening parameters for each band pass. 



